Mode mixing in Antiferromagnetically Correlated Double Quantum Wells 
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We examine the robustness of a recently predicted exchange- induced zero-field magnetic phase in 
semiconductor double quantum wells in which each well is spin-polarized and the polarization vec- 
tors are antiparallel. Magnetic instabilities are a general feature of Coulombic double quantum well 
systems at low densities. We argue that this antiferromagnetic phase is stabilized relative to ferro- 
magnetic ones by an effective superexchange interaction between the wells. Detailed self-consistent 
Hartree-Fock calculations using a point-contact model for the interaction show that the antiferro- 
magnetic phase survives intra-subband repulsion matrix elements neglected in earlier work in a large 
portion of the model's parameter space. We also examine the role of asymmetry due to biasing or to 
differences in the widths of the two quantum wells. The asymmetry creates a mode coupling between 
the intra- and inter-subband collective spin-density excitations (SDEs) which changes the Raman 
spectroscopy signature of the phase transition from a complete softening of the inter-subband SDE 
to a cusp as the density is tuned through the transition. This cusp may be detectable in inelastic 
light scattering experiments in samples of sufficient quality at low enough temperatures and densi- 
ties. 

PACS numbers: 73.20.Dx, 73.20.Mf, 75.70.Cn, 71.10.-w 



I. INTRODUCTION 

The subject of exchange-correlation-induced phase 
transitions has proven to be a rich field of research which 
has revealed many intriguing phenomena. Normal Fermi- 
liquid-state instabilities occur when the kinetic energy of 
the particles in a quantum system is of the same order 
or smaller than the inter-particle exchange and correla- 
tion energies. The instabilities lead to a variety of new 
electronic states. In Coulombic systems, this situation 
obtains at low densities and instabilities are more likely 
in systems of reduced dimensionality or, especially in two 
dimensions, in an applied magnetic field. Perhaps the 
best known examples of such new states occur in the in 
the fractional-quantum-Hall, strong-field regime of two- 
dimensional electronic systems .u 

The interest in exchange-correlation-induced phase 
transitions can be traced in part to the availability 
of high-quality semiconductor quantum wells, quantum 
wires, and superlattices. These artificial structures may 
be fabricated with remarkable precision and quality and 
possess electron densities that can be varied over a large 
range through a combination of modulation doping and 
judicious gating. The one- or two-dimensional nature of 
the resulting electronic systems, and the low densities 
realizable in devices of this kind, make them ideal for 
studies of interaction-driven physics. Another feature of 
these devices of importance to the current work is the ex- 



tra degree of freedom available when multiple layers are 
present, as in multiple quantum wells or superlattices. 
This new degree of freedom allows transitions to states 
with order not only in the intra-layer electronic degrees 
of freedom but also in diagonal or off diagonal inter-layer 
chargeoilj or spinBEJila observables. 

In particular, a great deal of attention has recently fo- 
cussed on the possibility of inter-layer spin ordering in 
wide single or double quantum wells. This attention is 
motivated by studies of quantum well structures in which 
the lowest two subbands are well separated in energy 
from the higher subbands and the density is sufficiently 
low so that only these two subbands are occupied. In 
these structures, an earlier theoretical calculation of the 
collective spin-density excitations (SDEs) in the absence 
of a magnetic field showed a complete softening of the 
inter-subband SDE in a range of densities around that at 
which the second subband begins to populate. til Subse- 
quent analysis indicated that this softening corresponded 
to a phase transition to a state in which each well in the 
double quantum well (or the electron gases on each side of 
a wide single quantum well) was spin polarized with the 
polarization vectors antiparallel; i.e., a transition-to an 
antiferromagnetic order in the well spin densities.cJ Fur- 
ther work predicted that a similar transition to a canted 
antiferromagnetic phase siould occur in the presence of 
a magnetic field at v = 2, til and this-jtransition seems to 
have been observed experimentally.113 However, the pre- 
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dieted transition in the zero-field case has not yet been 
observed. Ea 

Several possibilities exist which may explain the ab- 
sence of the zero-field antiferromagnetic phase in .these 
experiments. The calculations predicting this phaseta are 
based on a mean-field treatment of the interacting sys- 
tem, which is known to overestimate the densities and 
temperatures at which such symmetry-breaking transi- 
tions occur. In addition, electronic scattering by dis- 
order or impurities generally has a detrimental effect 
on correlation-induced phases. Both of these difficulties 
are exacerbated by the low dimensionality of the dou- 
ble quantum well system considered. Thus, the zero- 
field antiferromagnetic phase may exist but may not have 
been observed due to measurements made at tempera- 
tures which are too large in samples of insufficiently high, 
quality. The fact that the v = 2 transition is-observedt3 
as predicted by the mean-field calculationsEZl does not 
contradict this point of view, since the magnetic field 
completely quenches the kinetic energy and makes the 
mean-field theory a controlled approximation to the in- 
teracting system. 

Alternative explanations for the absence of the zerp* 
field phase lie in the structure of the theory itself,t£l 
which was derived to explore the qualitative features 
of the antiferromagnetic phase without considering sev- 
eral confounding effects which may nonetheless be impor- 
tant. In particular, the previous calculationEj did not ac- 
count for interactions between electrons in the same sub- 
band, which should be of the same size or stronger than 
the interactions between electrons in different subbands 
which were included. This omission can become espe- 
cially important when the intersubband excitation soft- 
ens. These interactions could introduce ferromagnetic 
or charge-ordered phases into the model which are not 
probed by the current experiments. More seriously, the 
slight asymmetry present in any realistic double quantum 
well structure will couple the intra- and inter-subband 
SDEs, potentially preventing the latter from softening. 
As this softening was expected to-ite-a hallmark of the 
antiferromagnetic phase transition,oE3 it seemed reason- 
able to hypothesize that any asymmetry in the structure 
might suppress the antiferromagnetic phase entirely. 

In this paper, we address the question of the robust- 
ness of the zero-field antiferromagnetic phase in the pres- 
ence of intra-subband interactions and asymmetry in the 
double quantum well structure. First, we argue on gen- 
eral grounds that the antiferromagnetic phase is a direct 
consequence of the importance of the intra- well exchange 
interaction at low densities and is stabilized by the inter- 
well hopping, which leads to an effective superexchangc 
interaction. Thus, this phase should obtain in a suitably 
constructed heterostructure. Second, we extend thepSfilf- 
consistent Hartree-Fock calculations of earlier worktj to 
include both intra- and inter-subband matrix elements of 
the model interaction and the effects of an asymmetric 
double quantum well. These calculations demonstrate 
that, while the intra-subband interaction docs introduce 



ferromagnetic phases and asymmetry does reduce the re- 
gion of the phase diagram occupied by the antiferromag- 
netic phase, the antiferromagnetic phase does not dis- 
appear. By examining the collective mode spectrum in 
the asymmetric structure, we also find that the inter- 
subband SDE does not soften due to the mode coupling 
between the intra- and inter-subband SDEs. The antifer- 
romagnetic transition nevertheless occurs as a result of 
the collapse of the mire-subband SDE which, through the 
mode-coupling, has a strongly antiferromagnetic charac- 
ter. 

The outline of this paper is as follows. In the next Sec- 
tion, we employ a simple model for two weakly coupled 
two-dimensional electron gases to examine the energetics 
of the antiferromagnetic transition. Sec. [II contains the 
formalism for the extended self-consistent Hartree-Fock 
theory used in the remainder of the paper. This for- 
malism includes asymmetry and all matrix elements of 
the interaction and is used to compute both the ground 
state and collective mode properties in what follows. In 
Sec. IV, we present the results of our computations for 



the matrix elements, phase diagram, and collective modes 
in this model and discuss their implications. Finally, in 
Sec. [y|, we summarize our results and conclude. 



II. ORIGIN AND STABILITY OF THE 
ANTIFERROMAGNETIC PHASE 

In this Section, we examine a simple model of a double 
quantum well in order to extract the basic physics un- 
derlying the zero-field antiferromagnetic phase. To that 
end, consider two two-dimensional electron gases sepa- 
rated by a barrier constructed so that the interaction 
between them is negligible. This would be the case if the 
two-dimensional layers were widely separated or the bar- 
rier were very high and if the dielectric constant of the 
barrier were large. Suppose the electrons in each layer 
move freely except for a Hubbard-like point-contact in- 
teraction V(v) = Vo<5(rL-jas used in textbook treatments 
of itinerant magnetism.Ej 

For equal charge density in the two layers, the Hartree- 
Fock energy -Ehf of the two-layer system can be derived 
following Refs. |,|JlO| and may be expressed in the form 
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where n is the total electronic density in both layers, Nq 
is the single-spin, two-dimensional density of states, and 
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is the relative spin polarization in layer i = 1,2 with 
partial spin-dependent densities rii a . In this equation, 
the first term represents the contribution to the energy 
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from combined kinetic and exchange effects, and the 
last term is the Hartree contribution. As is clear from 
Eq. (0), when N Vq > 1, it is energetically favorable for 
both layers to acquire a spontaneous spin polarization 
| m 1 1 = 1 771 2 1 = 1, while for NqVq < 1, the layer remain 
unpolarizcd. This result is .-simply the Stoner criterion 
for 2D itinerant magnetism.lia 

This description of spontaneous spin-polarization in 
two-dimensional electron systems is unrealistic both in 
the use of a point-contact interaction and in the use of 
the Hartree-Fock approximation. Our objective in this 
section is to obtain a qualitative understanding of the 
influence of weak electronic tunneling between the lay- 
ers on a double-layer system when isolated single layer 
systems are close to their ferromagnetic instabilities. We 
postpone a realistic discussion of the system parameters 
for which the physics we address in this paper is likely 
to realized to Sec. [v| For a sufficiently strong intra-layer 
repulsion, then, the exchange interaction forces both lay- 
ers to spin polarize, but the relative orientation of the 
polarizations is unknown. For simplicity, let us restrict 
our attention to two possibilities for the relative orien- 
tation: parallel (ferromagnetic) alignment or antiparallel 
(antiferromagnetic) alignment. With the spin unpolar- 
ized (paramagnetic) phase, the three possible phases for 
the two-layer system are shown in Fig. [l]. 

In the absence of any inter-layer coupling, the ferro- 
magnetic and antiferromagnetic phases are degenerate, 
but this situation changes if we introduce a small am- 
plitude for hopping between the layers. Writing Ciko- for 
the annihilation operator of an electron in layer i = 1, 2 
with two-dimensional wave vector k and spin projection 
a, this inter-layer hopping is governed by the Hamilto- 
nian 

#L = —J- [ ct lka C 2k CT 
ktr 

where Ao is the splitting between the symmetric and an- 
tisymmetric single-particle eigenstates. The leading or- 
der change in the ground state energy due to H± can 
be calculated using linear response theory, and is pro- 
portional to the transverse pseudospin susceptibility of 
double-layer systems defined in Ref. For the present 
model with no inter-layer interactions and contact intra- 
layer interactions, the time-dependent Hartree-Fock ap- 
proximation for the susceptibility gives 



h.c. 



(3) 
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~A 



(paramagnetic) 



4 -JV Aq (ferromagnetic) 
— |Aq/Vq (antiferromagnetic) 



(4) 



The fact that the results for paramagnetic and ferromag- 
netic states are independent of the interaction strength 
is a special property of the present model related to both 
the absence of interlayer interactions and the wave vector 
independence of the exchange self-energy. 

We see that, in all cases, introducing inter-layer hop- 
ping reduces the ground state energy, as one might ex- 
pect when the confinement of the electrons to the layers 



Paramagnetic 
Phase 



Ferromagnetic 
Phase 




Antiferromagnetic 
Phase 



Left 



Right 



Well Well 

FIG. 1. Population of spin and well states for the phases of 
a double quantum well considered in Sec. [n]. In the paramag- 
netic phase, there is no spin polarization in either well. In the 
magnetic phases, each well is completely spin-polarized with 
the polarization vectors parallel in the ferromagnetic phase 
and antiparallel in the antiferromagnetic phase. 
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is weakened. Comparing the energies of the ferromag- 
netic and antiferromagnetic phases, we also observe that 
the degeneracy between these phases is broken by the 
hopping term. Specifically, the antiferromagnetic phase 
is found to be more stable than the ferromagnetic phase 
if NqVq < 2, implying that the inter-layer hopping opens 
up a region of antiferromagnetic order between the para- 
magnetic and ferromagnetic phases. The mechanism for 
this stabilization can be deduced from the form of SE/A 
to be a superexchange interaction; that is, an electron is 
able to hop from one layer to the other and back at the 
cost of a Hubbard energy in the intermediate state, yield- 
ing an energy savings of Aq/2Vq. A similar mechanism is 
blocked by the Pauli exclusion principle in the ferromag- 
netic phase , since the hopping Hamiltonian preserves 
spin [Eq. (§]. 

We argue that this mechanism favoring an antiferro- 
magnetic arrangement of the ordered moments in the two 
layers will be dominant in most circumstances. The cal- 
culations described below indicate that an effective su- 
perexchange interaction between the wells should stabi- 
lize an antiferromagnetic phase for moderate interaction 
strengths. This explanation of the zero-field antiferro- 
magnetic phase suggests why the presence of asymme- 
try and additional interaction matrix elements may not 
eliminate this phase: both the two-dimensional ferro- 
magnetism within each well and the superexchange in- 
teraction between the wells should be fairly insensitive 
to these perturbations. In the remainder of the paper, 
we perform a more detailed self-consistent Hartree-Fock 
calculation to support this statement and to explore the 
consequences of these perturbations. 



III. FORMALISM 

To accomplish the goal of investigating matrix element 
and asymmetry effects, we employ an extension of the 
point-contact model described in Ref. |l5|. In the orig- 
inal calculation, the full Coulomb interaction between 
the electrons in the double quantum well was approxi- 
mated by a delta function in real space with only the 
matrix elements of this interaction between the lowest 
two subbands being kept and the remainder being set 
to zero.E-3 The use of a delta-function or point-contact 
interaction is based primarily on a desire to create a 
simple, solvable model which mimics the qualitative fea- 
tures of the fully interacting system. A quantitative the- 
ory would require that the interaction be made realistic, 
and also that the interactions be treated more accurately 
than in the Hartree-Fock approximation. In practise this 
would require quantum Monte-Carlo calculations of some 
type, which would involve an enormous amount of effort 
and would not be able to address the excitation spec- 
trum which provides the experimental signature for the 
state we are proposing. We therefore maintain the point- 
contact form of the interaction here. However, there is 




W 
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FIG. 2. Diagram of the asymmetric double quantum well 
considered in this paper indicating the well width w, barrier 
width b, confining potential U, and potential asymmetry dU. 
The values of the parameters used here and. m, previous studies 
of magnetic instabilities of these systemsEJ'Ej model a typical 
GaAs/AbGai-a, As double quantum well with w = 140 A, 
b = 30 A, and U = 220 meV. We consider both symmetric 
(dU = 0) and weakly asymmetric (dU = 0.5 meV) double 
quantum wells in what follows. 

no reason to set matrix elements of this interaction other 
than those between the lowest subbands to zero, as in 
previous work. Indeed, we shall see below that these 
other matrix elements are of the same order as the inter- 
subband ones. Thus, we shall include all the matrix el- 
ements between the lowest two subbands in our calcula- 
tions. 

In addition to the issue of intra-subband repulsion, we 
would also like to study the effects of quantum well asym- 
metry on the phase diagram and collective modes of this 
system. This asymmetry arises in real quantum wells 
through alloy fluctuations across the profile of the well 
or fluctuations in the well thickness, and we model it by 
allowing one of the wells to be deeper than the other, 
as illustrated in Fig. ||. Although we have assumed the 
effective asymmetry to enter through the well depth, the 
results of our calculations for systems with asymmetric 
well widths should be qualitatively similar. 

We therefore consider a three-dimensional electron gas 
confined along the z-direction by a potential Vconf{z) of 
the type shown in Fig. |^ which interacts through a three 
dimensional point-contact potential ^(R) = VS (R). In 
the absence of the interaction V^R), the electronic eigen- 
states are given by the solution of the time-independent 
Schrodinger equation 



2m* dz 2 



Vconf(z) 



(5) 



where m* is the effective mass, which is assumed to be 
constant throughout the heterostructure. Normalizing 
these eigenf unctions by 
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I dz\Uz)\ 2 = 1, (6) 
we write the electron annihilation operator ip a (R) as 

W(R) = -J= J2 e ik - r £n(^Kk CT , (7) 

where R = (r, z) = [x, y, z), k = (k x , k y ), A is the trans- 
verse area of the sample, and c„ko- annihilates a quasi- 
particle in subband n, of transverse wave vector k, and 
with spin projection a (these conventions will be used 
throughout this paper). 

Defining a composite subband and spin index a — 
(n a ,a a ) with summation over repeated indices implied, 
the Hamiltonian for this system is 

H = X! £ak C lk C ak 
k 

~^~2A ^ ^ ad ' bc C ak+q C bk'-q C ck' 
kk'q 

Here, the quasiparticle energy 

h 2 k 2 



c dk- 



(8) 



£ak 



2m* 



(9) 



is measured with respect to the chemical potential \x and 
the matrix elements of the interaction arc 



V ab , cd = 6 aa „ b 6 ac „ d V J dz& a (z)( nb (z)£ e (z)£ nd (z). (10) 

We treat this Hamiltonian within self-consistent 
Hartree-Fock theory allowing for the possibility of phases 
with broken symmetry in subband and spin indices, but 
imposing translational invariance within each layer. The 
electronic Green's function in the interacting system can 
therefore be written 



G a b{k n ) = - 







(T T [c k(r)4(0)]), (11) 



where k n = (k, iu> n ), f3 — l/T (h = kg = 1 throughout 
this paper), and the rest of the notation is standard.EiJ 
This Green's function is determined self-consistently 
from the self-energy in the Hartree-Fock approximation, 



[Va 



b,dc 



Vac 



Mb 



-iu?mO- 



Gcd(k'J, (12) 



and the Dyson equation, 

[fan - £ ak) S ab - ^ab] G bc (k n ) = Sac, 

under the constraint of constant sheet density N s , 



A 



Gaa (^m) j 



(13) 



(14) 



which determines the chemical potential. 

We solve these equations in the following way. Given 
a double quantum well structure defined by the effective 
mass to*, sheet density N s , and the structural parame- 
ters shown in Fig. |^, we compute the eigenfunctions and 
eigenenergies by solving Eq. (^|) with the normalization 
condition Eq. (^) . These eigenfunctions are used to com- 
pute the matrix elements of the interaction via Eq. ([To]) in 
terms of a single parameter V. The resulting matrix ele- 
ments and eigenenergies are employed to solve Eqs. (|TT|)- 
(|l4|) self-consistently at T = following the procedure 
outlined in Sec. II. B of Ref. [l^ and including all matrix 
elements of the interaction between the lowest two sub- 
bands. This procedure yields the ground state properties 
of the system as a function of the sheet density N S) the 
interaction parameter V, and the structural parameters 
of the double quantum well [Fig. 0]. 

To illuminate the properties of the interacting system 
further, we also compute the generalized density-density 
response function, which is defined by the analytic con- 
tinuation of 



n^(Q,w„ 



-iQR 



V jo 
(T T [^(R,r)^(O,0)]) 



(15) 



to real frequencies. In this expression, V is the system 
volume and the generalized density operator is 



^(R) = $>J( R )<C'^'( R ) 



(16) 



with ipa (R) given by Eq. (0) and = (l,a x ,a y ,a z ) are 
the Pauli matrices. This response function is computed 
from the non-interacting response function in subband 
and spin space within the conserving approximation de- 
scribed in Sec. V.A of Ref. [l5| but with the inclusion of 
all interaction matrix elements. In addition to the infor- 
mation this response function reveals about the excita- 
tions of the interacting system, its imaginary part is pro- 
portional to the intensity obsepjed in resonant inelastic 
light scattering measurementsJ23'Ej allowing us to make 
contact with experiment. This is particularly relevant 
here, because-searches for the antiferromagnetic phase 
in both finiteEj and zerolia magnetic field have employed 
this technique. 



IV. RESULTS 

In th is Section, we apply the formalism described in 
Sec. Ill to compute the ground state phase diagram and 
collective modes in a typical GaAs/ALjGai^As dou- 
ble quantum well structure which is expected to exhibit 
the zero-field antiferromagnetic instability. The structure 
has a well width of 140 A, a barrier width of 30 A, and 
a well depth of 220 meV [cf. Fig. ||], and an electronic 
effective mass to* = 0.067 m e . For the moment, we leave 
the asymmetry unspecified. 
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A. Matrix Elements 
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As a first step in obtaining the phase diagram for 
this structure, we must solve the time-independent 
Schrodinger equaution [Eq. (||)] for the lowest two eigen- 
functions £n(z) and eigenenergies e n at a fixed value 
of the asymmetry parameter dU [cf. Fig. || and then 
compute the matrix elements of the interaction through 
Eq. (|Io|). Solving Eq. (||) is straightforward and yields 
a splitting of the non-interacting eigenstates of Ag AS = 
£2 — £1 = 2.25 meV for dU — 0. As dU is increased, 
this splitting increases to a maximum of 18.5 meV at 
dU = 9.4 meV; for larger dU, the lowest eigenfunction 
is localized in one well. Since the structures examined 
experi ment ally have subband splittings on the order of 
1 meV,E3ll3 we restrict our attention to small dU values. 

The dependence of the matrix elements on the asym- 
metry parameter is somewhat more interesting and mer- 
its a brief discussion. Since our model interaction is 
a delta function in real space and we have chosen the 
wave functions to be real, the matrix elements defined 
by Eq. ( |Io|) are invariant under permutation of the in- 
dices. Thus, there are only five independent matrix ele- 
ments, Vu,n, V 2 2a2, Vn.22, Vn,i2, and V 22 ,2i, which are 
displayed in Fig. as a function of the asymmetry param- 
eter dU. In a symmetric double quantum well (dU = 0), 
the 1^11,12 and V22,2i matrix elements vanish by symme- 
try, but the remaining inter-subband (^11,22) and intra- 
subband (Vn^i and V22, 22) matrix elements are equal 
to within 5 %. We remark that for rnpre realistic inter- 
action models Vn.,22, which is roughlyO proportional to 
the difference of intra-layer and inter-layer interactions, is 
weaker than Vii,n and V22,22 which are roughly propor- 
tional to the sum. The latter matrix elements would be 
equal if the electrons were localized to the wells; the fact 
that they are nearly so indicates that the wave function 
overlap between the wells is small. Additionally, this cal- 
culation provides direct evidence that the neglect, of the 
intra-subband repulsion employed in earlier workEj is not 
generally justified for these double quantum well struc- 
tures. We shall see, however, that their inclusion in the 
calculation changes the qualitative picture only slightly. 

As the asymmetry is increased from zero, we discern 
several features. We see that the intra-subband matrix 
elements and V-J2.22 are approximately equal and 

increase with increasing dU to saturate at a value about 
twice the dU = one. The approximate equality of these 
diagonal matrix elements follows from the normalization 
condition, Eq. (^J), imposed on the two eigenfunctions. 
The increase in these matrix elements with dU, on the 
other hand, can be attributed to the increasing confine- 
ment of the wave functions of the two subbands to op- 
posite wells, similar to what occurs in the presence of 
an applied electric fieldB'Q Thus, at large dU, the two 
subband wave functions are almost completely localized 
in opposite wells, enhancing the magnitude of the diag- 
onal matrix elements. For the same reason, the inter- 
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FIG. 3. Dependence of the interaction matrix elements 
Vij.ki/V [Eq. (llQ)] on the double quantum well asymmetry 
parameter dU [cf. Fig. [|. Shown are the intra-subband ma- 
trix elements Vii,ii (dot-dashed line) and ^22,22 (dotted line), 
the inter-subband matrix element Vi.1,22 (solid line), and the 
asymmetry-induced matrix elements Vn,i2 (long dashed line) 
and V22,2i (short-dashed line). For the point-contact interac- 
tion employed in this paper, the matrix elements are invariant 
under permutation of the indices, so these matrix elements 
span the entire set. 

subband matrix element Vwpi decreases with increasing 
asymmetry: as the wave functions from different sub- 
bands are increasingly localized, their overlap, and hence 
Vii,22) decreases to zero. Since the zero-field antiferro- 
magnetic transition depends on this matrix element, it 
is clear that large asymmetry is inimical to this phase. 
Based on the small values of the observed splitting of-the. 
lowest two subbands in the experimental samples JiSEij 
however, we expect the actual samples to be in a regime 
in which Vn t 22 is still non-negligible. Hence, the antifer- 
romagnetic phase is not immediately excluded. 

Finally, we note that the mixing terms Vii,i2 and Vii,i\ 
have opposite signs and increase in magnitude with dU 
to a maximum around dU = 1 meV. The wave func- 
tion for the lowest (n = 1) subband has no nodes and 
we chose it to be positive. Orthogonality requires that 
the n = 2 wave function have a node and we choose it 
to be negative in the well where |£i(z)| 2 is largest; with 
this convention Vn 12 is negative and V22,2i is positive, 
as observed in Fig. |^. In addition, since these matrix ele- 
ments must vanish both in the symmetric [dU = 0) limit 
and when dU is large and the subband wave functions 
are localized in different wells, the maximum seen in this 
figure is also expected. These results suggest that the 
mode coupling between intra- and inter-subband excita- 
tions induced by these matrix elements will be maximal 
around dU = 1 meV. 

Taken together, the behavior of the interaction ma- 
trix elements presented in Fig. ^ indicate that the zero- 
field antiferromagnetic phase will probably not be stable 
against large asymmetry in the quantum wells. However, 
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the current experimental samples have subband split- 
tings more consistent with small asymmetry, and there- 
fore these samples may be of high enough quality to ob- 
serve this phase, at least in principle. To examine this 
situation further, we shall compute the phase diagram 
and collective modes for two choices of dll in the weak 
asymmetry regime: dll — (the symmetric case) and 
dll = 0.5 meV (the asymmetric case). The following 
subsections discuss the results of these calculations. 



B. Phase Diagram 

As described in the preceding subsection, the structure 
of the double quantum well yields the eigenenergies and 
eigenf unctions of the non-interacting system which are 
then used to compute the interaction matrix elements 
up to an overall factor V [Eq. (10)]. To be consistent 
with earlier work,E_j we choose to parameterize the inter- 
action strength by the magnitude of the inter-subband 
repulsion Vn,22 = V12 rather than by V, but all the ma- 
trix elements are uniquely determined by either param- 
eter. With the structure and interaction strength fixed, 
the only other parameter in our model is the sheet den- 
sity N s . Given these parameters, Eqs. (0)-(|l3) can be 
solved at T = to yield the interacting ground state of 
the system. The resulting phase diagrams in terms of 
the dimensionless interaction strength ./V0V12 and sheet 
density N s /2NoAg AS are presented in Fig. |^ for the dou- 
ble quantum well structure of Fig. 0. In these figures, 
No is the single-spin, two-dimensional density of states, 
and Ag AS = 62 — 61 is the subband splitting in the non- 
interacting system. 

Consider the symmetric (dU = 0) case first [Fig. §](a)]. 
At first glance, this phase diagram is very similar to 
the one obtained from the neglect of intra-subband 
rcpulsionjlj even though this repulsion is included in 
our calculations. At low density, the paramagnetic 
phase with one spin-degenerate subband occupied (la- 
beled Ni in the figure) is stable, while at higher densities 
(N S /2N A^ AS > 1) but weak interaction (N Vi 2 < 1/2) 
the paramagnetic phase with two spin-degenerate sub- 
bands occupied (labeled N2) is stable. For larger in- 
teraction strengths and intermediate densities, we see a 
broad region in which the antiferromagnetic (AF) phase 
obtains. Thus, the inclusion of intra-subband repulsion 
does not eliminate the zero-field antiferromagnetic phase 
from the phase diagram of the model, despite the fact 
that the intra- and inter-subband repulsions are of the 
same order [Fig. 

These matrix elements do have an effect on the phase 
diagram, however. At higher densities (N s /2NqAc iAS 
larger than approximately 2) and strong interaction 
(JV0V12 > 1/2), two different ferromagnetic phases ap- 
pear. In these phases, each quantum well is spin po- 
larized, and the polarization vectors are parallel. They 
differ in the arrangement and filling of the interacting 
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FIG. 4. Phase diagram of (a) a symmetric and (b) an asym- 
metric double quantum well in terms of the sheet density 
N a and the interaction parameter V12 for the point-contact 
model described in Sec. [oi|. The points are the numeri- 
cally computed boundaries between the various phases of this 
model which are labeled as follows: Ni and iV2 are param- 
agnetic phases with one and two spin-degenerate subbands 
occupied, respectively; FMi are ferromagnetic phases with i 
non-degenerate subbands occupied; and AF is the antiferro- 
magnetic phase. See Fig. |l| for a real-space depiction of these 
phases and Fig. ^| for the structural parameters of the dou- 
ble quantum well. In the figure, iVo is the single-spin, 2D 
electronic density of states, and Ag AS is the splitting between 
the lowest two states in the non-interacting system. Note that 
asymmetry reduces the extent of the antiferromagnetic phase 
in (b) but does not destroy it altogether. 
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bands, as indicated by the labels in the figure: the phase 
FMi corresponds to an interacting band structure with i 
spin-split subbands occupied. The presence of the ferro- 
magnetic phases is expected when intra-subband repul- 
sion in present; what is surprising is that the ferromag- 
netic phases do not exclude the presence of the antifcr- 
romagnetic phase. As argued in Sec. || based on a weak- 
coupling model, the antiferromagnetic phase is stabilized 
by a superexchange interaction for intermediate inter- 
action strengths. Similar behavior is seen in Fig. |J(a) 
based on our strong-coupling computations and presum- 
ably originates from the same mechanism. Note that, 
in the limit of vanishing hopping between the wells, V\2 
may be associated with Vq/2 in the model of Sec. |l im- 
plying that the phase boundary NqV\2 — 1/2 in Fig. |J(a) 
is nothing but the Stoner criterion for the related Hub- 
bard model. The fact that the model of Sec. [H] does not 
have a closer correspondence to Fig. ||](a) suggests that 
the relative magnitudes of the intra- and inter-subband 
interaction matrix elements, which are all equal in the 
model of Sec. H are important for determining whether 
ferro- or antiferromagnetic phases obtain in a specific re- 
gion of the phase diagram. 

With the introduction of weak asymmetry [dU = 
0.5 meV), the qualitative features of the phase diagram 
do not change, as seen in Fig. ^(b). As before, we find 
paramagnetic phases at small interaction strengths, the 
antiferromagnetic phase at larger interaction strengths 
and intermediate densities, and ferromagnetic phases at 
larger interaction strengths and higher densities. The 
identification of these phases precisely matches those 
in the symmetric case of Fig. ||(a), although the po- 
sition of the phase boundaries have shifted somewhat. 
An apparently new feature occurs at low density and 
large interaction strength, where an FM\ ferromagnetic 
phase has replaced the paramagnetic N\ phase. How- 
ever, this phase also occurs in the symmetric structure 
when N V\2 > 1, but is cut-off in Fig. ||(a). It corre- 
sponds to the usual ferromagnetic instability in a single, 
spin-degenerate band_syhich obtains when the interaction 
is sufficiently strong.Ej 

Although weak asymmetry clearly does not destroy the 
antiferromagnetic phase, it does have observable conse- 
quences. The most noticeable effect of the asymmetry on 
the ground state is that the spin polarizations in the mag- 
netic phases, which can be obtained from the expectation 
value of the density operator, Eq. (|li]), are no longer of 
equal magnitude in the wells. This is an obvious con- 
sequence of an asymmetric structure which nonetheless 
does not disturb the identification of ferro- and antifer- 
romagnetic phases, since one can determine whether the 
spin polarizations are parallel or antiparallel without re- 
ferring to their magnitudes. 



C. Collective Modes 

The effects of asymmetry on the antiferromagnetic 
phase cannot be fully appreciated based on the ground- 
state properties alone, but must be augmented by an 
examination of the excitation spectrum of the system. 
We focus on the spin-density excitations in what follows, 
since they are the excitations most strongly coupled to 
the magnetic nature of the ground state and cap .also 
be probed experimentally by Raman scattering. For 
these calculations, we compute the spin-spin response 
function (^i = 3 in Eq. (Ila)) for the appropriate ground 
state as discussed in Sec. pl| and identify the collective 
modes by peaks in the imaginary part of this response 
function. Since this procedure is used in the Raman 
scattering measurements, our results have direct impli- 
cations for experiment and we shall discuss them in this 
context. For concreteness, we fix the interaction strength 
and sweep the sheet density through the second-order 
transition from the N± to the AF phase in these calcula- 
tions [cf. Fig. f§. 

As an introduction to the general phenomenology of 
spin-density excitations in double quantum wells, con- 
sider the symmetric (dU = 0) case first. By appropriately 
arranging the light scattering geometry, Raman scatter- 
ing can selectively probe the intersubband spin-density 
excitations,E3E3 which, in our approximation, have the 
form shown in Fig. |^(a). In addition to a continuum of 
inter-subband particle-hole excitations, there is a collec- 
tive spin-density excitation (SDE) with a finite energy 
at q = which disperses with increasing q toward the 
particle-hole continuum. The magnitude of the q = 
SDE energy is reduced from the subband splitting A$as 
by vertex corrections appearing iu-the response function 
due the the exchange interaction.Ocj In addition, there 
is an intra-subband SDE which has a linear-in-q disper- 
sion in our model, shown by the dashed line in Fig. g|(a). 
If the well is symmetric, this mode will not appear in 
Raman spectra taken in a scattering geometry meant to 
observe inter-subband excitations. A symmetric system 
with identical quantum wells is invariant under inversion 
about the mid-point between the wells so that all states 
can be classified by a parity quantum number. Inter- 
subband excitations, which are odd, and intrasubband 
excitations, which are even, do not interact and can cross 
as seen in Fig. ||(a). 

As the density increases in our model, the exchange- 
induced reduction in the q = inter-subband SDE energy 
increases until the mode softens entirely, as illustrated 
in Fig. |^(b). This complete softening was seen initially 
in time-dependent, local density approximation calcula- 
tions of the SDE spectrum in these systems and was the. 
first evidence for the zero-field antiferromagnetic phase.O 
An analysis of the real-space spin response identified 
this softening as an antiferromagnetic transition of the 
well spin polarizations.li3 As the density is increased past 
the antiferromagnetic transition, the inter-subband SDE 
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FIG. 5. Collective spin-density excitations in a symmetric 
double quantum well of the type shown in Fig. | with dU = 0. 
(a) Collective mode spectrum in terms of excitation energy E 
and intra-well wave vector q showing the inter- (solid line) 
and intra- (dashed line) collective spin-density excitations as 
well as the continuum of inter-subband particle-hole excita- 
tions (shaded area). Although not apparent in the figure, the 
inter-subband collective mode does disperse with q. (b) q = 
inter-subband spin density excitation (SDE) energy E as the 
sheet density N 3 is tuned through the Ni to AF transition at 
fixed interaction strength N0V12 = 0.55 [cf. Fig. ^(a)]. Note 
that the SDE (solid line) softens completely at the transition 
point and that the collective mode appearing on the antiferro- 
magnetic side (A r s /2A r oAg AS > 0.919) is the amplitude mode 
of the antiferromagnetic order parameter (dashed line). The 
dot indicates the point in parameter space presented in (a). 
In this figure, iVo is the single-spin, 2D electronic density of 
states, Ag AS is the energy separation of the lowest two sub- 
bands in a non-interacting double quantum well, Asas is the 
splitting in the interacting system, and q\ — 2m*Ag AS /?i 2 . 
The intra-band particle-hole continuum is not indicated. 



turns into the collective mode associated with amplitude, 
fluctuations of the antiferromagnetic order parameter .113 
Experimentally, then, one expects to see a complete soft- 
ening of the inter-subband SDE and the recovery of this 
amplitude mode as the density is tuned through the tran- 
sition. 

The presence of asymmetry in the double quantum 
well complicates this picture somewhat. The SDE spec- 
trum for our asymmetric (dU = 0.5 meV) double quan- 
tum well illustrates these complications and is shown in 
Fig. ^(a). Most noticeably, the asymmetry mixes the 
intra- and inter-subband excitations, so that even in scat- 
tering geometries designed to measure only inter-subband 
response, both intra- and inter-subband excitations will 
appear.EJ This effect is seen through both an enlarged 
particle-hole continuum and the presence of a damped 
mode in the inter-subband spectrum corresponding to the 
intra-subband SDE. Furthermore, the asymmetry cou- 
ples the intra- and inter-subband SDEs themselves, lead- 
ing to an avoided crossing which may be seen by com- 
paring Fig. ||(a) to Fig. ||(a). 

This avoided crossing effectively prevents the inter- 
subband SDE from completely softening on entering the 
antiferromagnetic phase, as shown in Fig. ^|(b). In this 
figure, ones sees that, as the density is tuned toward 
the antiferromagnetic transition, the q = inter-subband 
SDE energy decreases to a finite value. Further increase 
of the density into the antiferromagnetic phase reverses 
this trend, and the energy of the amplitude mode of the 
antiferromagnetic order parameter increases, resulting in 
a cusp. 

Despite the fact that the inter-subband SDE does not 
completely soften, the antiferromagnetic phase appears, 
demonstrating that this softening is a sufficient but not 
necessary signature of the zero-field phase transition. We 
expect some mode to soften at this transition, however, 
and it turns out that it is the miro-subband SDE which 
collapses. Specifically, the mode coupling between inter- 
and intra-subband SDEs pushes the latter mode down 
in energy at all q, effectively reducing the group velocity 
of this mode. Approaching the antiferromagnetic transi- 
tion, the q — > group velocity of the intra-subband SDE 
decreases until it vanishes at the transition point. The 
resulting phase has the character of the real-space spin 
density profile of the intra-subband mode, which direct 
calculation reveals to be antiferromagnetic. The charac- 
ter of the intra-subband excitation at small q is therefore 
very similar to the inter-subband excitation in the sym- 
metric well due to the mode coupling between the two 
excitations. Thus, its softening can lead to an antiferro- 
magnetic transition without inconsistency. 



V. CONCLUSION 



In this paper, we have examined the origin and ro- 
bustness of a zero-field antiferromagnetic phase in dou- 
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FIG. 6. Collective spin-density excitations in a weakly 
asymmetric double quantum well of the type shown in Fig. 
with dU — 0.5 meV. The organization of and notation in 
the figure are the same as Fig. H with the interaction strength 
fixed at A r oVi2 = 0.55 and the antiferromagnetic transition oc- 
curring at N B /2NoAg AS — 0.698. As seen in (a), the asymme- 
try induces a coupling between the inter- and intra-subband 
spin-density excitations which results in the avoided cross- 
ing of inter- and intra-subband spin-density dispersions and 
a mixing of intra-subband particle-hole excitations with the 
inter-subband continuum. In (b), one sees that this mode 
coupling also prevents the inter-subband SDE from softening 
when going from the paramagnetic to the antiferromagnetic 
phase. 



ble quantum wells predicted in previous calculations 
Based on a simple model, we determined that such mag- 
netic phases are a direct consequence of the magnetism 
expected at low densities in two-dimensional systems 
where interaction effects dominate. In particular, mag- 
netic phases with either ferro- or antiferromagnetic spin 
polarizations are possible, and a superexchange interac- 
tion between the wells leads to a region at intermediate 
densities and interaction strengths where the antiferro- 
magnetic phase is preferred. 

By performing a detailed self-consistent Hartree-Fock 
calculation on a model of these double quantum well sys- 
tems, r*e addressed two features left out of preceding 
work:lia intra-subband repulsion and asymmetry of the 
heterostructure. Our results support and extend those 
of Ref. [l^. Specifically, we found that the antiferro- 
magnetic phase is stable in a large region of the model 
phase diagram despite the presence of intra-subband re- 
pulsion that is as strong as the inter-subband repulsion 
which drives the antiferromagnetic instability. The intra- 
subband repulsion does, however, introduce ferromag- 
netic phases, producing a rich phase diagram. Note, how- 
ever, that no charge ordering phases of the type discussed 
m Refs. g-O were observed, even though our formalism 
did allow for that possibility. In addition, both ferro- 
and antiferromagnetic phases persist in the presence of 
asymmetry in the quantum well structure; indeed, the 
phase diagram is qualitatively unaffected by its intro- 
duction. Asymmetry does have a strong influence on the 
collective mode spectrum, though, and induces a mode- 
coupling between inter- and intra-subband spin-density 
excitations which prevents the latter from completely 
softening at the antiferromagnetic phase transition. In- 
stead, the inter-subband spin-density excitation exhibits 
a cusp at the transition while the q — > group velocity 
of the intra-subband excitation vanishes. The coupling 
between these two modes nonetheless lends an antiferro- 
magnetic character to the intra-subband excitation in the 
asymmetric system and enables the collapse of its group 
velocity to yield the antiferromagnetic phase. 

Taken together, these results strongly indicate that, 
if an exchange-correlation-induced ferromagnetic transi- 
tion occurs in a single-layer system at sufficiently low 
density, then the antiferromagnetic phase in a double- 
layer system should also occur. However, the issue 
of whether or not the ferromagnetic phase obtains in 
a single two-dimensional layer is^-nat settled. Vari- 
ational Monte Carlo calculationsc3'CJ show the pres- 
ence of a ferromagnetic phase between the paramag- 
netic and Wigner crystal uhases, but Green's function 
Monte Carlo computationscJ find no intervening ferro- 
magnetism. More receal 1 -ijijurnerical work based on the 
Monte Carlo techniquet^'EirEI once again favors the exis- 
tence of a ferromagnetic transition. Other calculations 
using density- fi i ri rtional theory in the local-spin-density 
approximationi3iij also support the existence of a fer- 
romagnetic transition at sufficiently low density. Thus, 
although a definitive demonstration of ferromagnetism in 
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a single, two-dimensional layer is lacking, a large body of 
evidence exists which firmly supports this hypothesis. 

The precise conditions under which these phases would 
be observable are somewhat harder to elucidate based 
on the mean-field theory presented in this paper. The 
principle difficulty is that mean-field calculations will 
tend to overestimate the temperatures and densities at 
which exchange-correlation-induced phase transitions oc- 
cur. The source of this difficulty lies in the neglect of fluc- 
tuations in the theory, which play an important role in 
the low-dimensional structures considered. If one goes to 
densities and temperatures deep below the critical values 
of these parameters, we expect that the mean- field thepsy 
will give an accurate qualitative picture of the phase£3 

In sum, work on the single layer system indicates that 
its ferromagnetic state is unlikely to occur for electron- 
gas density parameters smaller than r s w 10, a much 
lower density than would be indicated by the Hartree- 
Fock approximation, for which the transition to the fer- 
romagnetic state occurs at r s w 2. The present work 
suggests that the two-layer antiferromagnetic phase, as 
well as two-layer ferromagnetic phases, are likely to be 
present in double-layer systems when the density per 
layer approaches the low value at which the single-layer 
ferromagnetic instability occurs. For the GaAs systems 
studied experimentally it therefore seems unlikely that 
the antiferromagnetic state will occur for densities per 
layer substantially larger than 10 10 cm -2 . However, it 
is exceedingly difficult to estimate the transition density 
theoretically and one must rely on experiment. 

Currently, a single experimental publication regarding 
a search for the,antiferromagnetic phase has appeared in 
the literature,E3 and the results are equivocal. The au- 
thors of this study report inelastic light scattering mea- 
surements of the long-wavelength inter-subband collec- 
tive spin-density excitations as a function of density in a 
double quantum well which was expected, on the basis of 
the original theoretical work, to show the antiferromag- 
netic instability.Ea Instead of completely .softening at a 
finite density, as predicted in earlier work," il 15 l the inter- 
subband SDE shows no dr.ainatic structure down to the 
lowest densities measuredJia These results could be ac- 
counted for in at least two different ways. The most likely 
explanation is that the electron density of the sample, of 
order 10 11 /cm 2 ,E3 is above the critical density for the 
antiferromagnetic transition. Alternatively, the calcula- 
tions in this paper demonstrate that slight asymmetry 
in the double quantum well will prevent a complete soft- 
ening of the inter-subband SDE and yield a cusp as a 
function of density. Since the energy of the amplitude 
mode in the antiferromagnetic phase is similar to that of 
the inter-subband SDE away from the critical density [cf. 
Fig. |6|(b)] , measurements at a closely spaced grid of den- 
sities may be required to detect this cusp. In addition, 
the cusp may be broadened by impurity or fluctuation 
effects which are beyond our mean-field theory, further 
increasing the difficulty of detecting the transition. Thus, 
the current experimental results cannot exclude the ex- 



istence of the zero-field antiferromagnetic state, and its 
robustness as demonstrated by the calculations in this 
paper leaves us optimistic that such a unique exchange- 
driven phase can occur in nature. 
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